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ABSTRACT 

We present the results of a systematic investigation of spectral evolution in the Z source GX 
349+2, using data obtained during 1998 with Proportional Counter Array (PCA) on-board the 
RXTE satellite. The source traced a extended normal branch (NB) and flaring branch (FB) 
in the color-color diagram (CD) and hardness-intensity diagram (HID) during these observa- 
tions. The spectra at different positions of Z-track were best fitted by a model consisting of 
a disk blackbody and a comptonized spectrum. A broad (Gaussian) iron line at ^ 6.7 keV is 
also required to improve the fit. The spectral parameters showed a systematic and significant 
variation with the position along the Z-track. The evolution in spectral parameters is discussed 
in the view of increasing mass accretion rate scenario, proposed to explain the motion of Z 
sources in the CD and HID. 

Key words: accretion, accretion disc - binaries: close - stars: individual: GX 349+2 - stars: 
neutron - X-rays: stars - X-rays: general. 



1 INTRODUCTION 

Low Mass X-ray binaries (LMXBs) containing neutron stars 
are broadly divided into two classes: Z and atoll sources 
(Hasinger & van der Klis 1989). This classification is based on the 
shape of the path traced by these sources in the X-ray color-color 
diagram (CD) or hardness-intensity diagram (HID). Z-sources trace 
out a Z-type pattern in CD and HID, and have luminosities close to 
Eddington limit (L^dd)- A Z-track consists of three parts, horizon- 
tal branch (HB), normal branch (NB) and flaring branch (FB). Atoll 
sources trace out a fragmented shape in CD. They have luminosi- 
ties 0.01-0. 1 Lg dd and probably, have weaker magneti c field s than 
Z-sources |vm der Klis 1995; Psaltis. Lamb & Mille r]|l995h . It is 
widely believed that in Z-sources the inferred mass accretion rate 
increases along the Z-track from HB to FB (Hasinger et al,. I99j]) . 
Similarly, in atoll sources inferred mass accretion rate increases 
along the fragmented path from the island state to upper part of 
the banana branch. 

X-ray spectra of LMXBs can be described as the sum of soft 
and hard components. Two different approaches have been used 
to model these components. In the first case, the soft compo- 
nent is modeled by a single temperature blackbody correspond- 
ing to the hot surface of neutron star and a hard component 
corresponding t o comptonized emission from the inner disk or 
boundary layer iWhite et alJll986t Iwhite. Stella & Parmadligsl 
iDi Salvo et all200(ll200lh ~ In the second approach, the soft com- 
ponent is described by a multicolor disk blackbody from accretion 
disk and the hard component is modeled by comptonized emis- 
sion (Mitsuda et al. 1984, 1989; Di Salvo et al. 2002). Comparison 
of these models with observed data have shown that both the ap- 
proaches require an additional broad iron Ka line, arising from the 
ionized disk or from the accretion disk corona (ADC). Atoll sources 



exhibit two distinct spectral states; namely hard and soft. In hard 
state their spectrum extends upto energies 100 keV I Bar ret et alj 
2000). The spectrum of Z sources is much softer, with cutoff en- 
ergies of a few keV. In addition, a hard tail has been occasionally 
seen in the spectrum of the 5 Z sources (GX 5-1 ' Asai et alJ^94 ; 
GX 17+2: Di Salvo et al. 2000; GX 349+2: Di Salvo et alj 2001 : 
Sco X-1: t D'Amicoe tSll2001l; Cyg X-2: iFrontera et aull99g : 
iDi Salvo et al. 2002) . 

UV/optical observations <VrtileketalJll99Cl) and correlated 
spectral and timing study of Z sources suggest {van der Klis 199^ 
that accretion rate increases along the Z-track from HB to FB. Mod- 
eling of X-ray spectra of Z sources also indicate that their motion 
i n the CD is caused by v ariations in the inferred mass accretion rate 
dPsaltis^ Lamb & Miller. ■1995..). However, other possibilities also 
exist ivan der Kli J2OOOI) Thus, investigations of spectral evolution 
along the CD may provide an important clue to understand this sce- 
nario. The evolution of X-ray spectra of sever al Z sour ces have been 
investi gated using Bep poS AX ( GX 17+2: Di Salvo et al. 2000; GX 
349+2: lDi Salvo et ali200L Cyg X-2: .Di Salvo et al..2002,) . In GX 
349+2, no appreciable change was observed in spectral parameters 
as source moved from NB to FB. Only, seed photon and electron 
temperature were slightly higher in FB than NB. Also, in GX 17+2 
spectral parameters did not vary significantly from HB to NB. X- 
ray spec tra of Cyg X-2 show ed an interesting evolution along the 
Z-track jDi Salvo et all2002h . It was found that as the source Cyg 
X-2 made a transition from HB to NB the inner disk temperature 
increased and at the same time inner rim of accretion disk moved 
inward. Also, optical depth of comptonized component systemat- 
ically decreased and electron temperature increased as the source 
moved from HB to NB. The authors suggest that hardening of the 
soft blackbody component and softening of the hard comptonized 
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Soft Color 

Figure 1. Color-color diagram of GX 349+2 for the January and September- 
October 1998 data sets. Soft color is count-rate-ratio in 3.5-6.4 keV and 2.0- 
3.5 keV energy ranges and hard color is that in energy bands 6.4-9.7 keV 
and 9.7-16.0 keV. The solid line in this figure represents an average Z-track 
passing through normal points. Each points is derived from a time bin size 
of 256s. 



emission are probably responsible for HB to NB transition. 
In this paper, we present the most detailed investigation of 
spectral evolution in GX 349-1-2 using extensive data set ob- 
tained with RXTE. The s ource traced the most complete 
CD during our observations (Zhana. Strohmaver & Swank 1998; 
lAgraw al & Bhattacharva. 2003; O'Neill et al. 2 002), providing an 
opportunity to study the spectral evolution in detail. 
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Figure 2. Count-rate spectrum and folded model ( power-law) for the Crab 
is shown in panel 1 . The residuals in the unit of a are shown in panel 2 (a 
systematic error of 1 % is included) and panel 3 (without an inclusion of 
systematic error). The filled circle represent the observed count rate spec- 
trum. 



2 OBSERVATION 

We analyzed public archival data obtained with PCA instrument 
on-board RXTE satellite. The data were collected during January 
and September-October, 1998. The total observation time with 
good data from PCA is ~ 220 ks. The PCA consists of five iden- 
tical proportional counter units wit h total effective are a of ~ 7000 
cm^ in the energy range 2-60 keV iTahodaetalll996l) . 



3 COLOR-COLOR DIAGRAM 

The PCA standard-2 mode data, accumulated every 16s in the en- 
ergy range 2-60 keV were used to create X-ray color-color diagram 
(CD) and Hardness-Intensity diagram (HID). The soft color is de- 
fined as the ratio of count rates in 3.5-6.4 keV and 2.0-3.5 keV 
energy ranges and hard color is defined as that in 9.7-16.0 keV and 
6.4-9.7 keV ener gy ranges. The 'rank number' or 'Sz' parameteri- 
zation technique rD ieters & van der Kli J200(i and reference there 
in) was used to define the position of the source in CD. We se- 
lected the normal points in the CD in such a way that they form a 
smooth curve. The color-color points in the CD are projected onto 
this curve. The Sz parameter for each projected point was calcu- 
lated by measuring their distance from NB/FB vertex. NB/FB ver- 
tex {Sz = 2) and end-point of the FB {Sz = 3) were taken as ref- 
erence points and rest of the Z curve was normalized according to 
the length of FB. The CD obtained from January and September- 
October observations of this source is shown in figure 1 . The solid 
curve in figure 1 represents the approximate Z-track. 



4 SPECTRAL ANALYSIS 

The standard-2 mode PCA data were used to create the spectra in 
the energy range 2.5-25.0 keV. The new cmbrighh'le model file 
"pca_bkgd_cmbrightvle_e3v20020201" for epoch 3 was used to 
calculate the PCA background. We used the observation (from 
obsid 30133-01-02-00 to 30133-01-02-08) of Crab Nebula to de- 
rive the level of uncertainty in the PCA response matrix. Fitting 
the Crab spectrum (in the energy range 2.5-25.0 keV) to a power- 
law with photon-index (F) = 2. 185±0.005 provided a reduced 
(Xi5 — "K^ /d-o.f) of 911/51. But, an inclusion of 1% systematic 
error reduced the value of xi to 41/51. Hence, a systematic error 
of 1% was added to all the spectra to take into account the uncer- 
tainty in the PCA response matrix. In figure 2, the effect of PCA 
systematics on the Crab spectrum has been shown. 

We divided the NB into 5 intervals and FB into 8 intervals 
to make an extensive investigation of spectral evolution in the 
source GX 349-1-2 along the Z-track. The source and background 
spectra within each interval were added separately to form a cu- 
mulative source and background spectrum. The background sub- 
tracted energy spectrum in the 2-25 keV range were fitted us- 
ing a two-component model. The multi-temperature blackbody 
(diskbb in XSPEC, Mitsuda et al. 1984) or a single tempera- 
ture blackbody (bbody in XSPEC) was used to describe the soft 
component. To model the hard component, comptonization mod- 
els compst rSun yaev & T itarchuk 1980) and comptt I Titarchu^ 
1994) of XSPEC were used. We find that diskbb-l-comptt model 
provided a better fit {y^ — 133/48) to the spectrum at upper NB 
(NBl) compared to bbody-l-comptt model {xl = 158/48). The 
diskbb-l-compst model provides a poor fit to the spectrum for in- 
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Figure 3. Count-rate spectram with folded model ( 
diskbb+comptt+gauss+pl) for the interval NB2 is shown in panel 1. 
The residuals in the unit of a are shown in panel 2 (a systematic error of 1 
% is included) and panel 3 (without an inclusion of systematic error). The 
filled circles represent the observed count rate spectrum. 




Figure 4. Variation of the best-fit diskbb parameters as a function of posi- 
tion (in terms of Sz parameter) along the Z-track. The inner disk tempera- 
ture Tf,f f is given in keV and the inner disk radius Reff is given in km. 
P'dbb (in gives percentage contribution of diskbb flux to the total flux 
(in the energy range 2.5-25.0 keV). 



terval NBl (xl ~ 643/49 ). It is to be noted that, in our case the 
temperature of diskbb component is nearly coinciding with lower 
energy range of the detector and hence diskbb and bbody model 
cannot be discriminated clearly. The analysis of X-ray bursts from 
GX 17-1-2 has suggested that the soft component c omes from the 
optically thick accretion disk i lKuulkers et alj2002h . Therefore, we 
chose diskbb-l-comptt model to describe the 2.5-25.0 keV spectra 
at various parts of Z-track. A Gaussian emission line at ~ 6.7 keV is 
required to improve the fit to the spectra for all the positions on the 
Z-track. Addition of this line decreased xl from 133/48 to 30.7/45 
for the spectrum of NB 1 ( probability of chance improvement of 
fit is ~ 2.3 X lO"'^*). The count rate spectrum for the source GX 
349-1-2 (during interval NB2) in the energy range 2.5-25 keV and 
best-fit to it are shown in the panel 1 of figure 3. The residuals in 
the unit of sigma, with and without inclusion of systematic error 
are shown in panel 2 and 3 respectively. 



5 RESULTS AND DISCUSSION 

The values of best fit spectral parameters are given in Table 1 (for 
NB) and Table 2 (for FB). The diskbb model represents the soft 
emission expected from the optically thick accretion disk and pro- 
vides two important disk parameters, color temperature (Tcoi) of 
inner disk and inner disk radius (Rin)- The inner disk radius Rin 
is expressed by: 

\kmj ^/cos9 



In this expression K is the disk blackbody normalization, D is the 
distance to the source in kpc, and 9 is the inclination of the ac- 
cretion disk. Distance of the source GX 349+2 is uncertain and 
probably l ies in the range 5-10 kpc ICooke"andPqnm3E^ 
[christian a nd SwankI (19971: IvanParadiis & McClintockl ll994 ). 
Optical observation suggests that if GX 349-1-2 is a LMXB 
( mass of companion < 2Mq) then inclination should 
he 8 > 25° for th e neutron star of mass < 2Mq 
^Wachter & Mareon ■l99d). X-ray properties of Z-sources suggest 
that there are two differe nt classes of Z-source s jKuulkers et al] 
ri994tlKuulkers. van der Klis and Vauahan 1996) : Cyg-like ( Cyg 
X-2, GX 5-1, Sco X-1) and Sco-like sources (Sco X-1, GX 
349-1-2, GX 17-1-2). One of the b est studied Cyg-like source Cyg 
X-2 has an inclination 6 ~ 60° iOrosz and Kuu lkers"l999l). As, 
it has been proposed that Sco-like sourc es are being viewed 
at lower inclination than Cyg-like sources teuulkers et anil994 
iKuulkers. van der KUs and VaughaiJ Il996l) . for the source GX 
349-1-2 inclination should be 25° < 8 < 60°. Therefore we as- 
sume a mean inclination 9 — 45° and distance D = 7.5 kpc for GX 
349-1-2 to derive its inner disk radius. It is to be noted that a different 
choice of inclination and distance in the range specified above will 
scale the inner disk radius by a factor of 0.59 to 1.58. However, as 
the inner disk radius for all the positions on the Z-track are scaled 
by the same amount, an overall trend of the change in the inner disk 
radius will not be affected. 

In the multicolor disk model the local emissi on is described 
by blackbody spectrum iShakura & Sunvaev Il973h . However, it is 
expected that in the inner region of the disk, opacity due to elec- 
tron scattering dominates over that due to absorption processes. In 
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this case the local spectrum can be significantly modified by comp- 
tonization and can be approximated by, 

/. = jzB4T,,i) (2) 

where, / — is the spectral hardening factor , Tcoi 

is the color temperature of inner disk, T^f / is the effective tem- 
perature of inner dis k, v is the frecjuency of radiation, Bi,{Tcoi) is 
the Planck function (Ebisawa et al. 1994). The effective radius and 
temperature of inner disk is given by, 

Reff=fR,n (3) 
Teff = ^ (4) 

By solving the vertical structure and ra diative transfer equa- 
tion of the disk, IShimura & Takahar3 il995h showed that for vis- 
cosity parameter a ~ 0.1, luminosity close to Eddington limit 
and mass of the compact object 1.4-10 Mq the soft emission from 
disk can be described by equation 2 with spectral hardening factor 
f=1.9±0.1. 

The variation in the parameters of soft component as a func- 
tion of Sz parameter is shown in figure 4. Our study revealed that 
Te/ / increased from ~ 1 keV to ~ 1.5 keV as source moved from 
upper to middle part of NB, i.e in the direction of increasing S,. 
Increase in Tef / can be interpreted in terms of inward motion of 
accretion disk with increasing accretion rate. In fact inner disk ra- 
dius Reff was found to decrease from upper NB ( ~ 10.5 km) to 
middle NB('~ 8.0 km). As source moved further from middle NB 
to upper FB, Te// decreased from ~ 1.5 keV to ~ 0.86 keV. At 
the same time Reff increased from ~ 12 km to ~ 30 km. This 
suggests that inner rim of accretion disk moves outward again with 
increasing mass accretion rate. 

As mentioned earlier, diskbb radiation is emitted by optically 
thick matter, which is in Keplerian flow. Since the radiation drag 
force exerted by radiation coming from the neutron star can remove 
significant angular momentum from the gas in Keplerian flow, ra- 
dial flow starts a few stellar radii above the neutron star surface. The 
inner edge of the accretion disk can be taken as the point at which 
disk-flow ends and radial flow starts. Note that the intensity of ra- 
diation will be attenuated exponentially (exp"'^'", where Tr is the 
radial optical depth) as it passes through disk. Therefore radiation 
drag cannot remove any angular momentum beyond a certain dis- 
tance (Rami) from the central compact obj ect and supersonic radial 
flow i s possible only within radius Rami J Miller. Lamb & PsaltisI 
Il998tl Miller & Lambll996h . This radius is obtained by setting Tr 
= 5 (since most of the radiation will be attenuated beyond this op- 
tical depth). Since radial flow probably begins at the inner edge 
of the accretion disk. Rami should be very close to Reff, i-e. 
Rami ~ Reff - Radial optical depth of accreting matter is given 
by equation 5 of Miller. Lamb & Psaltis ( 1998), which is 

rr = ^ (l-2M/r')-^ 
47rmp J 7r(r')7(r') r'h{r')v'^{r') 

In this expression, M is the mass accretion rate through the in- 
ner disk, ar is the Thomson scattering cross section, v^{r') is 
the radial velocity, h{r') is the half-thickness of disk, 7^ = [1 — 
(^f-)2j-i/2^^ ^ where is azimuthal veloc- 

ity of accretion flow. Note that in all the expressions of this section 
we have used units such that c = G = 1, where c is the speed of light 
and G is gravitational constant. By substituting Tr — 5 and solv- 



ing equation 1 for radius (see equation 6 of l Miller. Lamb & PsaltisI 
( 1998)) one can get 

^^ff-^+'(^^Y (^) (^) ^-(6) 

\0.01MeJ \W km J \0.1c J \10'2 J 

From equation 6, it is clear that inner disk radius will decrease 
with increasing accretion rate. This explains decrease of inner disk 
radius seen from upper NB to middle NB. However equation 6 can- 
not explain the increase of inner disk radius from lower NB to up- 
per part of FB. In this part of CD (from lower NB to upper FB) the 
source is supposed to be super-Eddington and hence gravitational 
pressure cannot support the radiation pressure. This will lead to an 
outflow in this part of the Z-track. It should be noted that equation 
5 and 6 do not take into account the presence of outflow. The pres- 
ence of outflow can remove the matter from radial inflow and can 
dump it into a hot central corona above the inner disk. This will 
lead to a reduction in the density of radial flow and hence, in radial 
optical depth. A decrease in the optical depth of radial flow will 
cause radiation to penetrate further inside the disk and hence the 
inner rim of the disk will shift outward (i.e, Ref / will increase). 
Also, if a fraction of matter from radial inflow is pumped in to a 
hot corona then the optical depth of this corona will increase. In 
fact optical depth of the corona is found to increase along the FB 
where most likely M > MecIcI- 

The contribution of the soft component (Pdbb) to the total flux 
increases from ~ 8% to 56% as the source moves from upper to 
lower NB. Again from lower part of NB to upper FB, the soft emis- 
sion decreased from ~56% to 17%. The variation in the soft com- 
ponent can be explained in terms of variations in the geometry of 
inner accretion disk with mass accretion rate. An inward motion of 
disk (from upper to middle NB) increased the area, emitting the soft 
component. As a result soft disk emission also increased. Similarly, 
an outward motion of disk (from lower NB to upper FB) caused a 
decrease in the soft disk flux. 

The comptonized spectrum is produced by upscattering of soft 
photons (described by the Wien spectrum at ~ 1 keV) by hot elec- 
trons in the central hot corona or boundary layer Our analysis 
showed that parameters of the comptonized spectrum evolved sig- 
nificantly as the source moved along the Z-track (see figure 5). The 
electron temperature Te of comptonizing plasma cloud increased 
from ~ 2.7 keV to 3.2 keV along the NB. At the same time optical 
depth r of the plasma cloud decreased from ~ 5.6 to 2.8. Prob- 
ably, increased heating of the central corona caused by increase 
in the incident flux reduces the density of corona and hence, de- 
crease r . A decrease in t leads to a reduction in relative energy gain 
of the comptonization process, given by comptonization parameter 
y = ^"^".j (where ruet? is the rest mass of the electron), from 
upper to lower part of NB. Other possibility is that a slow collapse 
of the central corona to a geometrically thin accretion disk, leaving 
low d ensity plasma cloud above the accretion disk tPi Salvo et alj 
l2002h may lead to a decrease in optical depth. In this case same 
amount of matter will be accreted on a smaller area (a strip around 
neutron star). This will cause an increase in seed photon temper- 
ature. In fact seed photon temperature Tw increased slowly from 
0.79± 0.01 keV to 0.90± 0.02 keV along the NB. However, it is not 
clear that why a central corona should collapse to a geometrically 
thin disk. As source moved further from NB/FB vertex to upper FB, 
Te decreased from ~ 3.3 keV to 2.6 keV and r increased from ~ 
3.2 to 6.4. In the FB, the radiation pressure is expected to be very 
high and dominant over gravitational pressure. Therefore, probably 
due to effect of radiation pressure, an outflow begins near the inner 
part of disk, creating a density enhancement in the hot corona above 
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Figure 5. Variation in tlie best-fit comptt parameters as a function of posi- 
tion (in terms of Sz parameter) along the Z-track. Electron temperature Te 
and seed photon temperature T^r are given in the unit of keV. 



the inner disk and hence causing optical depth to increase along the 
FB. The comptonization parameter y also increases along the FB 
due to increase in optical depth. 

The results of Di Salvo et J] <200ll) on the source GX 349-1-2 
from BeppoSAX data are different from our results. In their work, 
the whole Z-track was divided into two intervals; NB and FB. The 
spectra of FB and NB were best fitted with bbody-l-comptt-l-power- 
law-i-Gaussian-lines. As the source GX 349+2 moved from NB 
to FB temperature of the blackbody component slightly increased 
from 0.51±0.01 keV to 0.59±0.02 keV. At the same time opti- 
cal depth of the comptonizing cloud decreased from 11.7±0.4 to 
10.5±0.5 and electron temperature increased from 2.65±0.05 keV 
to 2.95±0.07 keV. The comptonization parameter y was ~ 2.84 
in the NB and ~ 2.54 in the FB. Our results show that diskbb 
component is harder in the NB (T^ff = 2.75 ± 0.01) than FB 
(Teff = 1.95 ± 0.01) and comptonized emission in the NB is soft 
iy ~ 0.38) compared to that in the FB {y ~ 0.63). There may be 
two possible reasons behind the differences seen between our 
results and those of Di Salvo et al. 1 2001). First, they have used 
a model which is different from that used in our work. Second, 
they have used different energy ranges to model the NB (0.1- 
200 keV) and FB (1.8-200 keV) spectra as LECS data were not 
available during the flaring state. 

We calculated the radius of seed photon emitting region (Rw) 

Rw ^ Sx lQ^D^[^/{Twf as given by 



using expression, Hw = 3 x WU^I Yf^/{Tw, 
lln'tZandet'ai]jl999l) . Here, D is the distance in kpc, Fhoi is the 
unabsorbed comptonized flux in ergs s~^cm^ , y is comptoniza- 
tion parameter and TV is seed photon temperature in keV. The 
Wien radius Rw was ~ 30-48 km in the NB and ^ 30-40 km 



in the FB. These values are considerably larger than that reported 
previously jPi Salvo etall2 00n. 

We found a strong broad iron line at energies 6.7 keV, with 
width ~ 0.7-1.4 keV. It is to be noted that iron line was stronger in 
the FB (eq. width = 154-295) than NB ( eq. width = 97-210). This 
behavior of GX 349-1-2 is just opposite to that exhibited by it during 
earlier observations (Di Salvo et al. 2001). We also note that in our 
case observed iron line is much stronger. 

Even though X-ray intensity decreases or remains almost con- 
stant along the NB, it is generally believed that accretion rate 
increases along the Z-track. This arguments mostly relies upon 
two observational facts. Firs t, UV and optical flux increases along 
the Z-track from HB to FB ( Vrtilek et al. 1990: ' van Paradiis et all 
[r99Qi) . UV and optical fluxes are t hought to be better in dicator of 
accretion rate than the X-ray flux iHasinseret ai]ll990h . Second, 
frequency of Horizontal-Branch-Oscillation (HBO) and kHz quasi- 
periodic-oscillation (QPO) increases from HB to NB ( van der Klli 
1200(1 and reference there in). In almost all the models, fre- 
que ncy of these two oscillations is dire ctly propo rtional to accretion 
rate jLamb et al. ll985trPsaltis et al. Il 999: Mill er. Lamb & PsaltisI 
Il998i: I Stella. Vietri & Morsink "l999). However, there are a few 
recent results which indica te that increasin g accretion rate scenario 
may not be the true story i Homan et al. 1l2002h and motion along 
Z-trac k is caused by variation in the inner disk radius ivan der KlisI 
l200d) . Our systematic investigation of spectral evolution showed 
that the motion of the source along the Z-track can be explained by 
variation in two emission components; soft disk emission and hard 
comptonized spectrum. The motion of the source along the NB (i.e 
softening of spectra) is caused by an increase in disk contribution 
and a decrease in degree of comptonization. Similarly, a decrease in 
the contribution of the soft component and an increase in Compton 
parameter y seems to be responsible for the motion of the source 
along the FB. As discussed earlier, the above changes in the soft 
component and in the comptonized spectrum are expected if accre- 
tion rate is increasing along the Z-track. Hence, our results seems 
to be in the favor of the scenario in which motion of the source 
along the Z-track is caused by variations in the accretion rate. 



6 CONCLUSION 

In this paper, we have carried out a detailed study of spectral evolu- 
tion in the source GX 349-1-2 using RXTE data. During our obser- 
vations the source is found both in the NB and FB. The continuum 
spectra at different positions of Z-track are well described by a disk 
blackbody and a comptonized component. We explain the evolu- 
tion of parameters of both soft and hard components in terms of 
variations in the mass accretion rate along the Z-track. This study 
suggests that a systematic increase in the strength of soft compo- 
nent and a decrease in comptonization parameter y are responsible 
for the motion of the source along the NB. Similarly, a decrease in 
soft disk emission and increase in the y parameter may cause the 
source to move along the FB. 
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X-ray spectral evolution of GX 349+2 



Table 1. Results of diskbb+comptt+gauss fit to the NB spectra of GX 349+2 in the 2.5-25.0 keV energy 
range. For each spectrum the range of Sz is given in the bracket. Errors on the best fit parameters are 

calculated hy using A = 1 (689?- confidence) 



Para. 


NBl 


NB2 


NB3 


NB4 


NB5 




(1.80±0.02) 


(1.84±0.02) 


(1.88±0.02) 


(1.92±0.02) 


(1.96±0.02) 


Tw (keV) 


0.79±().01 


0.825±0.008 


0.85±0.01 


0.88±0.01 


0.90±0.02 


Te (keV) 


2 69+°-°l 
^■"^-0.004 


2.822±0.004 


2.92±0.02 


S.OOiO.Ol 


3.20±0.01 


T 


5.62±0.02 


5.11±0.01 


3.33±0.03 


3,00±0.03 


2.80±0.03 


Rw (km) 


48.35±2.2 


44.07±1.41 


22.71±1.1 


29.97±1.35 


29.28±1.05 


y 


0.67±0.006 


0.57±0.003 


0.26±0.005 


0.22±0.005 


0.20±0.004 


Te//(keV) 


l.OSitO.Ol 


1.25±0.005 


1.54±0.003 


1.49±0.005 


1.45±0.005 


R^ff (km) 


10.47±0.58 


9.40±0.81 


8.04±0.21 


11.94±0.14 


12.67±0.26 


Fdbbi%) 


7.64 


14.7 


55.80 


56.10 


55.51 


El (keV) 


6.76±0.2 




6.80±0.2 


6.78±0.16 


6.80±0.14 


ai (keV) 


A 77+0.25 
^^■"-0.14 


0.61±0.34 


0.72±0.15 


n 7q+0-26 
"■'^-0.11 


0.681°:?? 


Eq. Width (eV) 


155 


97 


142 


184 


210 


xHd.o.f) 


30.7 (45) 


31.0(45) 


32.8 (45) 


31.0(45) 


36.0 (45) 



Table 2. Results of diskbb+comptt+gauss fit to the FB spectra of GX 349+2 in the 2.5-25.0 keV energy 
range. For each spectrum the range of Sz is given in the bracket. Errors on the best fit parameters are 
calculated by using Ax^ = 1 (68% confidence) 



Para. 


FBI 


FB2 


FB3 


FB4 


FB5 


FB6 


FB7 


FB8 




( 2.06 + 0.04) 


(2.14+0.04) 


(2.22 +0,04) 


(2.30 + 0.04) 


(2.38 + 0.04) 


(2.46 +0.04) 


(2.54 + 0.04) 


(2.64 + 0.06) 


Tyy (keV) 


0.91+0.01 


0.869 + 0.001 


0.88+0.01 


0.981+0.004 


1.024 +0.003 


1.03+0.01 


1.07+0.02 


+0.015 
-0.005 


Te (keV) 


3.27+0.05 


2.67+0.01 


2.66+0.01 


2.68+0.02 


2.69+0.01 


2.66+0.02 


2.635 +0.005 


2.61+0.01 


T 


3.18+0.03 


5.28+0.01 


5.23+0.02 


5.34+0.03 


5.63+0.02 


5.97+0.02 


6.08+0.02 


6.37+0.03 


(^) 


30.61 + 1.55 


38.21+ 1.75 


39.23+1.90 


32.95+1.55 


29.00+1.34 


29.69+1.45 


27.15+1.38 


31.35+1.61 


y 


0.26 +0.007 


.58 +0.005 


0.57 +0.005 


0.60+0.008 


0.67+0.005 


0.74 +0.008 


0.76+0.006 


0.83 +0.008 


T^jj (keV) 


1.40+ 0.01 


1.08 +0.005 


1.16+0.006 


0.99 +0.004 


1.01+0.003 


0.95 +0.004 


0.91+0.01 


0.86+0.01 




12.92+0.61 


15.65+0.25 


13.20+0.43 


18.85+0.41 


23.65+0.24 


26.62 +0.28 


27.90 +0.88 


29.59 +0.75 




47.80 


20.7 


18.8 


17.43 


26.90 


23.36 


22.39 


16.56 


El (keV) 


6.74+0.01 


6.73+0.06 


6.78+0.04 






6.94+Of 


7.02+0-3 


6.9+0.18 


a, (keV) 


0.76+0.05 


227 


0.89+0.06 


0.75+0.09 


0.87+0.1 


1.17+0.09 


1.21+0.2 


1.45+0.16 


Bq. Width (eV) 


207 


241 


154 


161 


222 


205 


295 




32.4 (45) 


40.1 (45) 


54.0(45) 


43.7 (45) 


35.5 (45) 


34.2 (45) 


50.9 (45) 


58.5 (45) 



